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Abstract 

An overview is given of the methods for treating complicated problems without 
small parameters, when the standard perturbation theory based on the existence of 
small parameters becomes useless. Such complicated problems are typical of quan¬ 
tum physics, many-body physics, physics of complex systems, and various aspects 
of applied physics and applied mathematics. A general approach for dealing with 
such problems has been developed, called Self-Similar Approximation Theory. A 
concise survey of the main ideas of this approach is presented, with the emphasis on 
the basic notion of group self-similarity. The techniques are illustrated by examples 
from quantum field theory. 
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1 Preliminaries 


It would not be an exaggeration to say that practically all interesting realistic problems 
cannot be solved exactly, so that one almost always has to resort to some methods of 
approximate solution. In obtaining a solution, it is highly desirable, before plunging to 
numerical calculations, to get an approximate analytical solution to the problem, which 
could help to understand the basic properties and specihc features of the considered case. 
It is just the possibility of deriving analytical presentations for approximate solutions that 
is our main concern in this paper. After deriving and studying such analytical presen¬ 
tations, nothing prohibits one to pass to a numerical procedure. Moreover, a prefactory 
analytical consideration can help in devising an approximate numerical algorithm, and the 
knowledge of the basic peculiarities of the problem under consideration can save plenty 
of computer time. 

A general approach for treating complicated real world problems has been developed, 
called Self-Similar Approximation Theory. In this paper, we aim at delineating the prin¬ 
cipal ideas of the approach and in presenting its several new developments. To clearly 
distinguish the pivotal aspects of our theory from the characteristic points of other known 
techniques, we feel it necessary to say several words on the latter. There are, roughly 
speaking, three common ways of obtaining approximate solutions: single step estimates, 
asymptotic perturbation theory, and methods of successive iteration. 


A. Single Step Estimates 

This kind of estimates is often related to minimizing or maximizing the corresponding 
part of an inequality. Probably, the most known and widely used such a tool is based 
on the Gibbs-Bogolubov inequalities, which are formulated as follows. Let A and B be 
Hermitian operators on a Hilbert space, for which the form 

F[A] = -TlnTre"^^ 


exists, where j3T = 1, with T being real. Dehne an average of A with respect to B as 


< A >B = 


Tie-^^A 

Tie-hB 
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and similarly, an average of B with respect to A. Then the Gibbs-Bogolubov inequalities 
are 

<A-B>a< F[A] - F[B] <<A-B>b ■ (1) 

This is the generalized presentation of the inequalities which are better known for the case 
when, instead of arbitrary Hermitian operators, one deals with Hamiltonians H and Hq. 
In that case, the right-hand side of Eq. (1) becomes the inequality for the free energies. 


F[H] < F[Ho] + <H-Ho>o, 


where < ... >o implies the averaging with respect to Flo. This inequality was derived by 
Gibbs [1] for classical statistics. Bogolubov [2] generalized it for quantum statistics and 
added the left-hand side of Eq. (1). When T —0, the free energy reduces to the ground- 
state energy. Then, as a particular case, one has the Peierls inequality. Introducing an 
effective Hamiltonian 

Heff = Hq < H — Ho >0 , 
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one may also write 


FIH] < F[H,„] . 

In the standard way, one chooses the approximating Hamiltonian Hq = de¬ 

pending on a set of parameters or functions, say on a trial spectrum Uq, so that Hq could 
model the considered system and would allow one to calculate the free energy F[Heff] for 
the effective Hamiltonian -He// = Heff{oJq). Then, one minimizes F[Heff\ with respect to 
the trial functions uiq given by the equation 

-F = 0 . ( 2 ) 

When < H — Hq >0 is positively dehned, an approximate minimization can be done with 
the help of equation 

< H — HQ{uq) >0= 0 . (3) 

The Gibbs-Bogolubov inequality is constantly used in various problems of statistical 
mechanics and condensed matter theory. Among thousands of examples, let us mention 
the self-consistent photon approximation [3]. In quantum mechanics, the minimization of 
an effective ground-state energy with respect to trial parameters incorporated in a trial 
wave function is usually named the Ritz variational method. Such methods are termed 
single step since they give just a single approximation, without hinting on how to obtain 
subsequent corrections. 

B. Asymptotic Perturbation Theory 

Contrary to the single step estimates, perturbation theory is a systematic procedure 
dehning a sequence of approximations. It may, of course, happen for some complicated 
problems that one technically is able to calculate only a few terms of the perturbation 
sequence, but perturbation theory is systematic in the sense of prescribing a general way 
for calculating perturbative terms of arbitrary order. There is, especially in physical 
literature, quite a mess in the usage of the terms ’’perturbation theory” as opposed to 
’’nonperturbative methods”. Therefore, to avoid in what follows linguistic confusion, it is 
necessary to concretize several principal points and definitions. 

The standard perturbation theory presupposes the existence of small parameters per¬ 
mitting one to present solutions in the form of asymptotic series [4]. Because of the 
latter, the standard perturbation theory may be named asymptotic perturbation theory. 
Usually, when talking about perturbation theory, one keeps in mind exactly the standard 
asymptotic perturbation theory. 

One can distinguish three main types of asymptotic series occurring in perturbation 
theory. Suppose the problem is in calculating a real function f{x) of a variable a; G X C M. 
The case of one function of one variable is taken just for the simplicity of notations and is 
not principal. In general, the function f{x) can depend on any number of other variables, 
but we separate and explicitly write down only one variable that is assumed to play the 
role of a small parameter. Perturbation theory with respect to a small parameter |a;| 1 

yileds a sequence of approximants <Pk{x), where /c = 0,1, 2,..., approximating the sought 
function f{x) in the vicinity of a; = 0, 

f{x) ~ (pk{x) (|a;| < 1) . (4) 
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The approximants (pk{x) have the structure of asymptotic series of one of the following 
types. 

(i) Expansion over small parameters: 


k 

‘Pkix) = OnX” . 
n=0 

Generally, this can be an expansion over one or several parameters, 
(ii) Expansion over asymptotic sequences: 

k 

^kiY) 'y ' • 

n=0 


Here, the sequence is asymptotic in the sense of Poincare, so that 


£n+l{x) 



0 ). 


In particular, en{x) can be e^-^x), with e{x) being a given function of x. 
(iii) Generalized asymptotic expansion: 


(5) 


( 6 ) 

(7) 


k 

^k{x) = Y an{x)en{x) . (8) 

n=0 

In the latter, the coefficients a„(a;) retain dependence on x in order to satisfy some addi¬ 
tional conditions, but so that the sequence {a„(a;)e„(a;)} be asymptotic. 

The general feature of all the cases above is that the difference 


Aipk{x) = ipk{x) - ipk-i{x) 


(9) 


forms an asymptotic sequence {A99fc(a;)}, such that 


AqPk+i{x) 

A(pk{x) 


{x 


0 ). 


The construction of generalized asymptotic expansions can often be rather elaborate, 
aiming at improving the accnracy of calculations. For example, in the Lindstedt-Poincare 
method [5], one expands over a small parameter the sought solution and the freqnency 
choosing the expansion coefficients so that to kill secular terms. In the Krylov-Bogolubov 
averaging technique [6], the generalized asymptotic expansion is constrncted above an 
initial approximation inclnding nonlinearity in order to model well the main properties of 
a nonlinear system. In the theory of anharmonic crystals [7,8], the expansion in powers of 
a small anharmonic parameter starts with a self-consistent phonon approximation partly 
taking acconnt of anharmonicity. Nevertheless, no matter how elaborate is an initial ap¬ 
proximation and how complicated is the strncture of the resulting generalized expansion, 
all abovementioned cases presuppose the existence of small parameters and are typical 
examples of asymptotic perturbation theory. 
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C. Methods of Successive Iteration 

For numerical iterative algorithms, the existence of small parameters is not required. 
When the considered equation can be presented in the form 

Af{x) = h{x) , 

in which A is an operator and h{x), a given function, then an iterative solution of this 
equation reads 

iPk+i{x) = (pk{x) + Bk [A(pk{x) - h{x)] , ( 10 ) 

where Bk are operators chosen so that to make the sequence {(pk{x)} convergent [9]. The 
iterative algorithm (10) leads to approximants (pk{x) that, in general, are not asymptotic 
series, even if a; —0. This is because an expansion of (pk{x) at x = 0, provided it exists, 
has the structure of a series 

Mx) = , 

n=0 

with coefficients labelled by two indices. Then the difference (9) is 

H (fln - X"" + an , 

n=0 

which is not of order x^ but contains lower powers of x, since ^ a^^. Hence, {A99fc(a;)} 
is not an asymptotic sequence. 

Though the method of successive iteration looks more general than perturbation the¬ 
ory, it has a weak point of being, in the majority of cases, purely numerical, without 
providing analytical formulas. And also, for very complicated problems a numerical pro¬ 
cedure can be too much time consuming or even unsolvable. 

2 Optimized Perturbation Theory 

A qualitatively different systematic approach was advanced [10] for treating the problems 
without small parameters and permitting one to obtain, at least in first orders, approxi¬ 
mate analytical solutions. This approach was motivated by and applied to many-particle 
systems with strong interactions [11-21]. Technically, the approach is based on the meth¬ 
ods of perturbation and iteration theories combined with optimal control theory. 

The pivotal idea of the approach [10] is the introduction of control functions whose role 
is to govern the convergence of approximation sequences. Generally, an approximation 
theory consists of three parts: an initial approximation, a calculational algorithm, and a 
sequence analysis. Control functions can be introduced in any of these parts. The main 
is that the resulting approximation sequence {Fk{x,Uk)} be convergent due to control 
functions Uk = Uk{x). Optimized approximants are 

fk{x) = Fk{x,Uk{x)) {k = 0,1,2,...) , (11) 

forming a convergent sequence {fk{x)}. In obtaining the approximations {Fk{x,Uk)}, it is 
not required to have any small parameters, but, if the techniques of perturbation theory 
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are employed, an expansion is made with respect to a formal parameter e that at the end 
is set to nnity 

k 

Fk{x, Uk) = ank{x)e'^ (e ^ 1) . 

n=0 

This expansion is not asymptotic, even when a; —0, since the coefficients Unk are labelled 
by two indices. In this way, the expressions Fk{x, Uk), althongh can be obtained by means 
of formal techniqnes of perturbation theory, have the structure typical of the terms of 
an iterative procedure. Thus, the approach is not an asymptotic perturbation theory. 
Moreover, the approximations Fk{x,Uk) may also be obtained by means of an iterative 
algorithm. Therefore, the optimized perturbation theory is principally different from the 
standard asymptotic perturbation theory, since the former does not require any small 
parameters and the structure of its approximate terms is not that of asymptotic series. 

Definition. Optimized Perturbation Theory is a systematic method dehning a se¬ 
quence of successive approximants, whose convergence is governed by control functions. 

How could we formulate general rules for hnding such control functions? Since the 
role of the latter is to govern convergence, let us write down the Cauchy criterion of 
uniform convergence for the sequence {Fk{x,Uk)}, with a; G X. The sequence is uniformly 
convergent on X if and only if for each given e there exists such that 


\Fk+p{x,Uk+p) - Fk{x,Uk)\ < E 


( 12 ) 


for all k > Nfp and p > 1. We also wish that convergence be as fast as possible. A 
general way of dehning control functions, providing for the considered system the required 
property, is formulated by the optimal control theory [22]. According to this theory, control 
functions are given by minimizing a cost functional. In our case, the cost functional 
garanteeing the fastest convergence can be constructed as the fastest-convergence cost 
functional 

OO 

F u X 'y ' [Ffi-\-p(^x, Ufi^pj Fjii^x, Ujfj] , (13) 

^ n=0 

in line with the Cauchy criterion (12). The fastest convergence is achieved by means of 
control functions, minimizing the cost functional (13), that is, from the condition 


6uk 


= 0 , 




> 0 . 


The variational derivatives are 


6Fu 

Suk 


{2Fk — Fk+p — Fk-p)Fl , 


(14) 


^ = 2{Fff + {2Fk - Fk+p - Fk.p)Ff , (15) 

where Fk = Fk{x, Uk), I < p < k, and 

duk' ^ dul ■ 
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The first of Eqs. (14) possesses two solutions yielding two possible optimization conditions, 
(i) Differential optimization condition 


— Fk{x,Uk) = t) (16) 

OUk 

gives an extremum of the cost functional, but it is not clear what this extremum is, since 
the sign of the second derivative 




{2Fk — Fk+p — Fk-p)Fl! 


is not dehned. To concretize the situation, we need to invoke some additional constraints. 
Consider a particular case of p = 1 and assume the validity of the fast convergence 
condition 

Fk{x,Uk) ^ Fk+i{x,Uk+i) . (17) 

Then the differential condition (16) is an approximate condition for the minimum of the 
cost functional, since 


6ul 


0 , 



6ul 


^ QiFkf > 0 . 


(ii) Difference optimization condition 

Fk+p — 2,Fk + Fk-p = 0 

clearly correspond to the minimum of the cost functional, as far as 




2(F',f > 0 . 


(18) 


However, this does not uniquely dehne control functions, since Eq. (18) contains three of 
them, Ufc+p, Uki and Uk-p- To resolve the problem, consider again the case of p = 1 and 
assume that the value Fk{x,Uk) weakly depends on the change of Uk by Uk+i, which can 
be formulated as the weak sensitivity condition 


Fk{x,Uk) ^ Fk{x,Uk+i) . (19) 

Then the difference condition (18) reduces to the equation 

Fk{x, Uk) - Fk-i{x, Uk) = 0 , (20) 

making it possible to define Uk = Uk{x). 

In this way, neither Eq. (16) nor Eq. (18) can serve as exact equations for minimizing 
the cost functional (13) and for uniquely dehning control functions. But the latter are 
unambiguously dehned and the cost functional (13) is approximately minimized by either 
optimization condition (16) or condition (20). These conditions are equivalent to each 
other. Both of them are approximate. Both are variational, following from the variation of 
a cost functional. Both invoke the notion of weak sensitivity with respect to the variation 


7 



of control functions. Both are designed for supplying fastest convergence of the sequence 
of optimized approximants by minimizing the fastest-convergence cost functional. Being 
completely equivalent by their meaning, the optimization conditions (16) and (20) dif¬ 
fer only by their form, being expressed either through a derivative or through a hnite 
difference. 

Solving one of the optimization equations, either (16) or (20), one gets control functions 
Uk = Uk{x) governing convergence of the optimized approximants (11). It may happen 
for some k, that an optimization equation does not have an exact solution. Then, since 
the optimization conditions themselves are approximate, it is admissible to look for an 
approximate solution for a control function. For example, if the latter becomes complex, 
one may take only its real part. Or, when neither the derivative (16) nor difference (20) 
are exactly zero, one may look for control functions minimizing one of the following forms: 

mm\Fk{x,Uk) - Fk-i{x,Uk)\ ■ 

Uk 

As is mentioned at the beginning of this section, control functions can be incorporated 
at any step of the theory. A straightforward way is to introduce them in the initial 
approximation. For instance, if one considers a problem described by a Hamiltonian H, 
one may take for the initial approximation a Hamiltonian Hq{u) containing a set of trial 
parameters u. Then the Hamiltonian of the problem can be presented as 

H = Ho{u) + e[H - Ho{u)] , 

with a formal parameter e —> 1. Accomplishing perturbative calculations with respect 
to the formal parameters e, one sets it to unity and obtains perturbative terms Fk{x,u). 
After this, one hnds control functions from an optimization condition and substitute them 
into Fk{x, Uk). Absolutely the same procedure can be realized if the considered problem is 
characterized not by a Hamiltonian but by a Lagrangian or an action. It is also possible 
to incorporate trial parameters into an initial approximation for a wave function or for 
Green functions and to derive the subsequent approximations by means of an iterative 
procedure. 

As is evident, there can be a variety of particular technical ways of introducing initial 
approximations and control functions. But all such variants are based on the same fun¬ 
damental idea of control functions governing convergence for the sequence of optimized 
approximants [10]. Several years after Ref. [10], there appeared a series of papers [23- 
28] advertizing the same idea of introducing control functions for rendering perturbation 
theory convergent. Nowdays the optimized perturbation theory is widely used for vari¬ 
ous problems, being employed under different guises and called by different names, such 
as modihed perturbation theory, renormalized perturbation theory, variational perturba¬ 
tion theory, controlled perturbation theory, self-consistent perturbation theory, oscillator- 
representation method, delta expansion, optimized expansion, nonperturbative expan¬ 
sion, and so on [23-32]. Many problems of quantum mechanics, statistical mechanics, 
condensed matter physics, and quantum held theory are successively treated by this opti¬ 
mized approach. Not to list numerous applications, let us cite a couple of recent reviews 
[33,34]. 

Despite a number of very successful applications, optimized perturbation theory does 
not provide answers to the following important questions: 


min Fk{x,Uk) 
Uk duk 



(i) How to improve the accuracy with a given number of approximate terms? 

(ii) How to realize a step-by-step control of the stability of the method, if no exact 
solutions for the problem are known? 

(hi) How to choose the best initial approximation, if several of them are admissible? 
To answer these questions it is necessary to look at the problem from the point of view 
of a more general approach, which is described in the next section. 

3 Self-Similar Approximation Theory 

An approach, more general than optimized perturbation theory, has been developed [35- 
57], being based on the techniques of the theory of dynamical systems and optimal control 
theory. The underlying idea of the approach is to consider the passage from one succes¬ 
sive approximation to another as the motion on the manifold of approximations. Then, 
the approximation order k = 0,1,2,... should play the role of discrete time, and each 
approximant should represent a point of a trajectory in the phase space of approximants. 

To correctly describe evolution, one has to construct a dynamical system. For this pur¬ 
pose, we again need to introduce control functions. Employing some variant of the method 
of successive approximations, we get the terms Fk{x,Uk)- The role of control functions 
Uk = Uk{x) is to make convergent the sequence {fk{x)} of the optimized approximants 

fk{x) = Fk{x,Uk{x)) . (21) 

Convergence implies the existence of the limit 

lim fk{x) = f*{x) . (22) 

/c—>oo 

Let all approximating functions fk{x), with k = 0,1, 2,... and a; G X, together with the 
limit (22) pertain to a complete space A, which we call approximation space and which 
plays the role of a phase space for the evolution of fk{x) with respect to the discrete time 

k. 

Introduce a function Xk{^) by the reonomic constraint 

Fo{x,Uk{x)) = (p , X = Xk((p) . (23) 

Changing the variables, we dehne 

= fkixki^)) ■ (24) 

The transformation ^ ^ ^ is an endomorphism of the phase space A, with a unitary 

element given by the equation 

yo{(p) = (p . (25) 

By dehnition (24), each fk{x) corresponds to yk{<p) and, conversely, to each yk{<p) we set 
in correspondence 

fkix) = ykiFo{x,Ukix))) . (26) 
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By this construction, the sequences {Vki,^)} and {fk{,x)} are bijective. The existence of 
the limit (22) implies the existence of the limit 

lim yk{(f) = , (27) 

/c—>oo 

so that 

r{x)=y*{Fo{xy{x)) , (28) 

where u*{x) = lim^^oo ^^(a;). For the map {?/*:(</?)}, the limit y*{(p) is a hxed point, when 

yk{y*{,v)) =y\v) ■ (29) 

The dependence of the hxed point y*{^) on the starting point (p is due to the reonomic 
constraint (23). 

A particular form of the endomorphism (24) depends on the choice of an initial ap¬ 
proximation Fq{x,u) and of control functions Uk{x). These are to be chosen so that to 
guarantee the fastest convergence of the sequence {yk{y^)}- Uniform convergence on A 
implies the validity of the Cauchy criterion 

\yk+p{p>) - yk{p>)\ < e (30) 

for all k > and p > 1. This suggests that the evolution in the map {pfc(<7’)} is to be 
such that to minimize the fastest-convergence cost functional 

2 OO 

= [yn+p{^) - yn{^)f ■ (31) 

The minimization is to be done with respect to any yn{p>), including |/o(<p) = ‘P- As is 
evident, the absolute minimum of the functional (31) is realized if and only if the initial 
point ip coincides with the hxed point |/*(<p), when 

J-, = 0, ip = y*{p>). (32) 

Proposition. For the fastest-convergence cost functional to be minimal it is necessary 
that the self-similarity relation 


yk+p{p>) = yk{yp{p>)) (33) 

be valid. 

Proof. The absolute minimum of the cost functional (31) is Fy = 0. This is realized 
if and only ii ip = p*(<p), when equation (33) becomes an identity. 

Clearly, the self-similarity relation (33) is a necessary but not sufficient condition 
for minimizing the cost functional (31). This relation describes the property of self¬ 
similarity in the general sense, which includes as a particular case the scaling self-similarity 
l/fc(A(p) = A"''pfc((p) taking place for the power-law functions Pfc(<7’) = with the powers 
such that ak+p = Ofc ■ Op. The scaling self-similarity is what one usually keeps in mind 
when mentioning this property. However, this type of self-similarity is just a particular 
trivial case of the relation (33). The latter for the endomorphism together with the 
unitary element (25), dehnes the semigroup properties 

yk+p yk' yp 1 yo i ^ (34) 
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because of which the relation (33) can be called the group self-similarity. In the the¬ 
ory of dynamical systems [58,59], the family of endomorphisms {yk\ k G Z+j, where 
Z_|_ = {0,1,2,...}, is termed a cascade, which is a dynamical system in discrete time. 
Since in our case this family of endomorphisms is formed of the sequence {yk{^)} of the 
approximants yk{<p), the corresponding cascade can be named the approximation cascade. 

To deal with discrete time is less convenient than with continuous time, when the 
latter is given on M+ = [0, cx)). Therefore, it is useful to embed the cascade into a flow. 
The embedding 

{yk\ k e C {y{t,.. .)\t e R+} (35) 

implies that the flow possesses the same semigroup property 

y{t + t',(p) =y{t,y{t',^)) (36) 

and the trajectory {y(t, (p)} of the flow passes trough all points of the cascade, 

yit, if) = ykip) {t = ke'L+) . (37) 


The flow embedding the approximation cascade is termed the approximation flow. 

The advantage of dealing with a flow is that differentiating the self-similarity relation 
(36) one comes to a differential Lie equation 

^^y{t,p) = v{y{t,p)) , (38) 


where viy) is a velocity held. Integrating the evolution equation (38) between yk-i = 
yk-i{p) and y^ = yl{p), where yl is an approximate hxed point, we get the evolution 
integral 


dy 


fVk 

'vk-i v{y) 




(39) 


with Tk being the effective approximation time required for reaching the quasihxed point 
yl. For short, we may call Tk the control time. Due to the relations (23) to (28), the 
integral (39) may be presented as 


rfk dip 

Ifk-i Vk{p) 


— 1 


(40) 


where fk = fk{x), ff = fl{x), with 


fki^) = y*kiFo{x,Mx)) , 


(41) 


and Vk{p) is the cascade velocity given by a discretization of the how velocity. The Euler 
discretization of the how velocity is dehned as 


Vk{p) = Vk{xk{p)) 


where 


with Uk 


Vk{x) = Fk{x,Uk) - Fk-i{x,Uk) + {uk - Uk-i) 


d 

duk 


Fk{x,Uk) , 


Uk{x) and k = 1,2,.... 


(42) 

(43) 
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Our aim is to find a fixed point y*{(p) of the map {yk{H^)}i which, by construction, 
corresponds to the sought function f{x). The hxed point for a flow is given by the zero 
velocity v{y*) = 0. For the approximation cascade, we have 


Vk{x) = 0 . 


(44) 


This condition is to be treated as an equation for control functions Uk{x). However, 
the cascade velocity (43) contains two control functions, Uk{x) and Uk-i{x), hence one 
equation (44) cannot dehne them both. Thus, we either have to invoke some additional 
constraints or can use an approximate equation for hnding control functions, for instance, 
by minimizing the absolute value of the cascade velocity 

|14(a;)| < \Fk{x,Uk) - Fk-i{x,Uk) \ + 


d 

{Uk - Ufc-i)7^— Fk{x, Uk) 
UUh 


(45) 


There are three ways of formulating hxed-point conditions. 

Fixed-point condition 1. Find ui{x) from an additional condition, say, from the 
differential optimization condition (16) or from the difference optimization condition (20). 
Then, all other Uk{x), with k >2, are given by Eq. (44), that is, by the equation 

d 

Fk{x,Uk) - Fk-i{x,Uk) + {uk - Ufc-i)7^— Fk{x,Uk) = 0 . (46) 

OUk 

The disadvantage of this way is that control functions for k = 1 and k >2 are defined by 
different conditions. Also, if for some k equation (46) has no solutions, then an ambiguity 
arises requiring some additional assumptions. 

Fixed-point condition 2. Minimizing the first term in the right-hand side of Eq. 
(45), one has 

Fk{x,Uk) - Fk-i{x,Uk) = ^ , (47) 

which coincides with the difference condition (20). Then the cascade velocity (43) is 

d 

Vk{x) = {uk - Uk-i)T \— Fk{x,Uk) . (48) 

UUh 


Since Eq. (47) gives Uk{x) starting from /c = 1, the function Uq{x) is left undehned. 
Hence, the velocity (48) is valid for k > 2. If Eq. (47) has no solution for some k, then 
one could hnd Uk{x) by minimizing \Fk — Fk-i\. But in that case, the cascade velocity is 
not given by Eq. (48). 

Fixed-point condition 3. Minimize the second term in the right-hand side of Eq. 
(45), which yields 


(Ufc — Uk-l) 


d 

duk 


Fk{x,Uk) = 0 . 


(49) 


This is to be understood as the equation 


d 

duk 


Fk{x,Uk) = 0 , 


'^k 7 ^ 1^k—l ) 


(50) 


provided a solution for Uk{x) exists, or as the equality 


— ^k—1 ’) 


d 

duk 


Fk{x, Uk^ 7 ^ 0 , 


(51) 
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when the differential condition (50) does not possess a solution for Uk{x). In all the cases 
under condition (49), the cascade velocity becomes 


Vk{x) = Fk{x,Uk) - Fk-i{x,Uk) , 


(52) 


where Uk = Uk{x) and k > 1. 

Comparing the possible hxed-point conditions (46), (47), and (49), we see that the lat¬ 
ter is more general and provides a unique unambiguous way for dehning control functions 
Uk{x) for all A; > 1. In all three cases, the function uo{x) is, generally, not dehned, so is 
the term fo{x) = Fq{x,Uo{x)). Therefore, the evolution integral (40) has to be considered 
starting with k = 2. The effective time in Eq. (40) can be treated as another control 
function. By its dehnition, is the minimal time required for reaching a quasihxed point 
from an approximant fk-i- In general, the minimal time should correspond to one 
step, which means that kxk should be close to one. Therefore, the effective approximation 
time can be evaluated as 

Tfc = ^ . (53) 

In this way, an approximate hxed point fk{x), representing a k-th order self-similar ap¬ 
proximant for the sought function /(x), is completely defined by the evolution integral 
(40). 

An important advantage of the self-similar approximation theory is the possibility 
of controlling the stability of the procedure, which can be done by invoking the ideas of 
dynamical theory. For this purpose, for a map {|/fc(v^)}, one may dehne the local multipliers 


The multiplier at a quasihxed point is given by 

Ki'p) = Mvlif)) ■ 


(54) 


(55) 


The images of these multipliers on the manifold X are obtained by means of the change 
of variables (23). The image of Eq. (54) reads 


mk{x) = fj,k{Fo{x,Uk{x))) , 

and that of Eq. (55) is 

ml{x) =pkifkix)) . 

A quasihxed point is stable provided that 


(56) 

(57) 


< 1 , \m.l(x)\ < 1 . (58) 

It is useful to consider uniform stability characterized by the maximal local multipliers 

ml = sup \ml{x)\ . (59) 

X 

Then a k-th order approximant is uniformly stable if 

yl<l, ml<l. (60) 


hi = sup \fili(p)\ , 
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Generally, because of the definition of the multipliers (55) to (57) through the same 
local multiplier (54) given on the phase space A, the maximal multipliers (59) coincide, 
fil = ml, so that it is sufficient to consider one of them. 

The stability analysis also makes it possible to answer the question ’’which of several 
admissible initial approximations should one prefer in calculating higher-order approx- 
imants?” The answer is straightforward: One has to prefer that initial approximation 
which garantees the best stability of the procedure [57]. 

In conclusion to this section, let us emphasize a principal aspect distinguishing our 
approach from different variants of the standard renormalization-group techniques [60-62]. 
In the latter, one tries to establish either an exact or an approximate relation between 
a function f{x) and its value f{\x) for the scaled physical variable. Such a relation 
describes the motion with respect to the scaling parameter. Contrary to this, in self-similar 
approximation theory, we do not scale physical variables. But the group self-similarity 
describes an evolution in the phase space of approximants, with the approximation order 
playing the role of time. 

4 Method of Fractal Transforms 

One of the basic ideas in self-similar approximation theory is the introduction of control 
functions which govern the evolution of an approximation dynamical system to be close 
to a hxed point. As was mentioned earlier, control functions can be introduced at any 
part of calculational procedure. For instance, this can be done at the step of choosing 
an initial approximation, which results in the sequence of optimized approximants. But 
this can also be accomplished in the last part of calculations, after deriving a sequence of 
perturbative terms. 

As is discussed in Section 1, employing the standard perturbation theory, one obtains 
approximations having the structure of asymptotic series. Such series are usually diver¬ 
gent and have no sense for hnite values of expansion parameters. There exist the so-called 
resummation methods ascribing hnite values to divergent series [63]. The most often used 
among such techniques are the Borel summation [63] and the construction of Bade ap¬ 
proximants [64], including the two-point [65] and multivalued [66,67] Bade approximants. 
These techniques have many known limitations. Thus, to get a good accuracy, they re¬ 
quire to invoke a number of perturbative terms which often are not available. And also, 
such techniques are, actually, numerical. 

In order to incorporate control functions into a given asymptotic series, one needs to 
resort to a transformation including some trial parameters. The transformation involved 
must decode the self-similarity property hidden in the given perturbative sequence. For 
power series, it looks natural to employ the power-law transformations [68-75]. Since 
power laws are typical of fractals [76,77] the power-law transformation can also be called 
the fractal transformation [75]. 

For a function f{x), the fractal transform is 


F{x,s) = x^f{x) , 

with a real s. The inverse transform is 

(61) 

f{x) = x~^F{x, s) . 

(62) 
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The fractal transform satisfies the scaling relation 


F{Xx,s) ^ F{x,s) 

f{\x) f{x) 

reminding us those typical of fractals. 

Assume that for a hnite function f{x), the standard perturbation theory in powers of 
X results in a set of approximations 

k 

^k{x) = ’ ( 63 ) 

n=0 

where >^o{x) = (pox'^^ 7 ^ 0 and a; —0. The powers are arranged in the ascending order 

7n < 7n+i (n = 0,1, 2,..., A; - 1) ; (64) 

the signs of 7 „ being arbitrary. The structure of the series (63) is that of an expansion 
( 6 ) over the asymptotic sequence en{x) = a;'*'". 

Since it is always more convenient to work with dimensionless quantities, we assume 
that the variable x is dimensionless. Dehne the dimensionless scale-invariant function 


gk{,x) 

_ P>k{x) 

~ ‘fo(x) ' 

( 65 ) 

Employing the notation 

_ 

= ) 

9^0 

0 

1 

III 

e 


we have 

k 


gk{.x) = 

-- ^ a^a;"" , 

( 66 ) 


n=0 


where 



go{x) = Qq - 

= 1 ttO = 0 ; 

( 67 ) 

and the powers are such that 




(n = l,2,...,k-l) . 

( 68 ) 

Introduce the fractal transform 



Fk{x,s 

) = x"gk(x) 

( 69 ) 

and the inverse one 



gkix) = 

x~^Fk(x, s) . 

( 70 ) 

For the series (66), this gives 



M 

= anx''^'^" . 

( 71 ) 


n=0 


Then we apply the same procedure of self-similar approximation theory, described in the 
previous section, to the sequence {Fk{x, s)}. The difference is that here the trial parameter 
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s is treated as a quasi-invariant, that is, it is kept fixed under the evolution of the map 
{Fk (x, s)} and it is transformed in a control function Sk{x) aposteriori, the latter being 
defined from convergence and boundary conditions. 

The reonomic constraint (23) now becomes 


Fo(x,s) = ip, x = x(ip,s). 


(72) 


With the form (71), from where Fo{x,s) = x^, this gives x{(p,s) 
phism (24) now is 


which results in 


yk{^,s) = Fk{x{^,s),s) , 
yk{^,s) = • 

n=0 


For s being a quasi-invariant, the cascade velocity (42) is 

Vki^, s) = Vki^, s) - Vk-ii^, s) . 
From Eqs. (74) and (75), one gets 

Vk{(p,s) = . 


The endomor- 

(73) 

(74) 

(75) 


This is to be substituted in the evolution integral 


rvk 


dy 


Vk{y,s) 


— '^k 1 


(76) 


which is similar to the integral (39), and where yl = yl{(p,s) is a quasifixed point, with 
yQ{{p,s) = ip and k>l. Note a slight difference between the integrals (39) and (76). In 
the present case, the evolution integral (76) is obtained by integrating the Lie equation 
(38) between two quasifixed points, y^-i and y^. After changing variables according to 
the constraint (72), the integral (76) reduces to 


rn 


dp 


’Fl-i 


— '^k 1 


where F^ = F^{x,s), k>l and 


= Vl(Fo(x,s),s) . 


(77) 


(78) 


For /c = 0, since y^^p, s) = p, one has Fq[x, s) = x®. 

Calculating the evolution integral (77), with the cascade velocity (75), results in the 
iterative equation 


(Ftf- = (Fuf + At , 


(79) 


in which 


Ak — OjkdkFk 


Ok = - 


s 

If we recall that, to return to the original physical quantities, we have to accomplish the 
inverse fractal transform (70), then we define 


gk{,x,s) = X "F^{x,s) , 


(80) 
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with go^XjS) = 1. For the transform (80), the iterative equation (79) reads 

gt = sfii + . (81) 

At this stage, we have to convert the trial parameter s into a control function Sk = 
Sk{x). Taking this into account, we come to a self-similar approximant 


9*k{x) = gk{x,Sk) . 

Equation (81) for the approximant (82) can be written as 


Sk 


9*kix) = Wgl-iix) 


i/Sk 


(82) 


(83) 


where 6k = —ak/sk and 


9o{x) = 1 . 


(84) 


The solution to Eq. (83), with the initial condition (84), gives gk{x), which dehnes the 
self-similar approximant 

fkix) = Mx)9*kix) (85) 


for the sought function f{x). The quantities Ak and 6k are expressed through control 
functions Tk and Sk, which actually means that Ak and 6k can be considered themselves 
as control functions. The latter are to be dehned from additional conditions, such as 
convergence and boundary conditions. 


5 Self-Similar Root Approximants 

When the behaviour of the sought function is known in the asymptotic vicinity of two 
boundaries of the domain X, the control parameters Ak and 6k can be found from the re¬ 
lated boundary conditions [71,73,74]. Here, we generalize this procedure to the case when 
the boundary asymptotic expansions contain arbitrary powers of x, including noninteger 
powers. 

Suppose, for concreteness, that the variable x is given on the real semiaxes X = M+ = 
[0, cx)). If this is not so, then it is always possible to resort to a change of variables reducing 
the domain of x to M+. As earlier, we keep in mind that the variable x as well as the 
sought function are normalized to dimensionless units, so that the considered function is 
presented in a scale-invariant form g{x). Assume that the asymptotic behaviour of g{x) 
in the vicinity of the left boundary, 

gix) ~ gkix) 

is given by the asymptotic series 

k 

9k{x) = Y. 

n=0 

where 

go{x) = ao = 1 , ao = 0 , (88) 


[X 


0 ), 


ar,X 


( 86 ) 


(87) 
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and the powers are arbitrary real numbers arranged in the ascending order 

0 < On < ctn+i (n = 1,2,..., A; - 1) . (89) 

And let us assume that the asymptotic behaviour of g{x) at the right boundary is also 
known, 

g{x) ~ Gk{x) {x^ oo) , (90) 

being presented by the asymptotic series 

Gk{x) = bnX^'^ , (91) 

n=l 

in which 

^ 0 , A ^ 0 , (92) 

and the powers are arranged in the descending order 

Pn>Pn+i {n = 1,2,... ,k - 1) . (93) 

Note that Pn can be of any sign. 

Iterating k times Eq. (83), and using the notation 

= ^ (p= 1,2, 1) , Ufc = ^ , (94) 

dp Ok 

we obtain the self-similar root approximant 

gl{x) = (... ((1 + A^x'^Y^ + A 2 x'^^Y + • • • + AkX^^^Y^ , (95) 

which can also be called a nested root or superroot. The control parameters Ap and Up, 
with p = 1,2,..., k, are to be dehned from the asymptotic coincidence of expressions (95) 
and (91) at the right boundary, that is from the asymptotic boundary condition 

g*k{x) ~ Gk{x) {x ^ oo) . (96) 

In this way, the crossover formula (95) extrapolates the sought function from the left 
boundary a; —0 to the whole interval [0, cx)). 

As is easy to observe, reexpanding Eq. (95) in small x ^ 0 does not reproduce 
the structure of gk{x) in the series (87). For this to hold, one, hrst, should require 
that ap = pa, and even then the expansion coefficients would not coincide with Op for 
p > 1. However, there is no need to demand that the asymptotic behaviour of g*{x) be 
identical with gk{x). Vice versa, such a restriction would essentially spoil the accuracy 
of the approximant for large x. It is important to bring to mind that the main aim 
of the self-similar approximation theory is to construct accurate expression uniformly 
approximating the sought function in the whole interval of the variable x G [0, cx)). The 
asymptotic expansion (87) at a; —0 are used only as constructing blocks. In general, 
it would be possible to modify the form of the self-similar approximant so that it could 
exactly reproduce gk{x) at x ^ 0. But this, from the point of view of an accurate uniform 
extrapolation, is neither necessary nor correct. 
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At the same time, by the asymptotic condition (96), the control parameters Ap and 
Up are dehned so that the asymptotic expansion of at a; —> cx) exactly coincide 

with the series (91). Such a construction, as is evident, can be reverted in the following 
sense. One could derive self-similar root approximants starting from the right asymptotic 
expansion (91) and htting the corresponding control parameters so that the expansion of 
the derived root approximants at a; 0 be coinciding with the series (87). Nevertheless, 
the construction of crossover formulas from the left to right seems more preferable because 
of the following. The region of validity for the expansion (87) is |a;| 1 and that for 

the expansion (91) is |a;| ^ 1, that is, the region of validity of the right expansion is 
essentially larger than that of the left expansion. This conclusion is conhrmed by a 
number of particular cases demonstrating that the radius of convergence of the series (87) 
is usually zero, while that of the series (91) is hnite. 

In order to dehne the control parameters Ap and Up from the asymptotic condition 
(96), one has to know how to present an asymptotic form of the root approximant (95) 
at X ^ oo. To the hrst glance, the procedure of obtaining an asymptotic, as a; —cx), 
expression from the superroot (95) looks ambiguous, since the powers Up are yet not 
known, hence it is not clear how to classify larger and smaller terms. In particular cases, 
really, there can be several ways, depending on the relation between the values of ap and 
Up, of classifying asymptotic terms of the superroot (95) at a; —> cx. This ambiguity can 
be overcome by requiring the uniqueness and generality of the procedure. 

Definition. The self-similar root approximant (95) is called to be uniquely dehned 
by the asymptotic condition (96) if and only if all control parameters can be uniquely 
determined from a general rule, whose form is invariant with respect to the values of ap 
and which is valid for arbitrary p = 1 , 2 ,.... 

Theorem. The self-similar root approximant (95) is uniquely dehned by the asymp¬ 
totic condition (96) if and only if the powers Up are given by the equations 

otpUp — Op+i = const , 

apHp apj^\ j3};—p -\- , ( 97 ) 

aknk = l3i {p=l,2,...,k-l). 

Proof. When x —> cx, it is convenient to introduce the small parameter e = x~^. In 
terms of the latter, the superroot (95) can be identically presented as 

gl{x) = (l + 5;. (i + (1 + ... 

+^ 25 “— (1 + B.e-rT ■ ■ T'T ’ ( 98 ) 

where 

1 A'^p 

Bi = — , Bp+i = {p=l,2,...,k-l) . 

/ii ^p+i 

To prove sufficiency, let us assume that Eqs. (97) hold true. Then, 

Clp+l Pk—p /3k—p+l ■ 
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Due to the descending order (93) of (3n, one has 

ctp+i — apUp = const > 0 . 

This unambiguously dehnes the classihcation of powers of e in the form (98). Expanding 
Eq. (98) in powers of e, we observe that the hrst k terms of the expansion coincide 
with all k terms of the series (91), that is the asymptotic condition (96) is satished. To 
prove necessity, we assume that the superroot (95) is uniquely dehned by the asymptotic 
condition (96). This implies, according to the above dehnition, that there exists a unique 
general expansion of the form (98) in powers of e, which is invariant with respect to ap 
and p. As is evident from expression (98), such a unique general expansion is possible if 
and only if ctp+i — apHp = const > 0. Then Eq. (98) can be unambiguously expanded in 
powers of e. This expansion is to be compared with the series (91) that can be written as 

G,{x) = bix”' il + ^ s*-* 6 + h 6 + ,,, + ^ i 

Comparing the hrst k terms of the expansion of Eq. (98) with the series (91), we obtain 
Eqs. (97). 

This theorem makes it possible to apply the same general rules for constructing various 
crossover formulas. Let us stress that the theorem is a new result that has not yet been 
published. 

Till now, we have considered the situation when the order k of the left expansion (87) 
coincides with that of the right expansion (91). How could we proceed if these orders 
were different? 

If the number of available terms from the left is less than that from the right, this 
is not as important, provided we know the law prescribing the values of which is 
usually known or can be easily guessed. This is because we, actually, do not need to 
have all coefficients a„, which are incorporated in the control parameters An, and these 
are determined through the coefficients bn and powers f3n of the right expansion (91). 
Therefore, when only are available, nevertheless, we may add to the left expansion the 
required number of terms up to the order of the right expansion. 

When the number of terms in the right expansion is less than that of the left expansion, 
the situation again is not dangerous. Say, the left expansion gk{x) is of order k, while the 
right expansion Gm{x) is of order m < k. In that case, we iterate Eq. (83) till gl_n^^i{x) 
and at the next step, we set gk-m{x). Iterating in this way Eq. (83) m times, we come to 
the self-similar root approximant 

+Ak-m+2X'^^-^+^T^-^+^ + ... + AkX^^)""^ , (99) 

where the notation (94) is used. The superroot (99) for the case m = k returns to the 
form (95), since go{x) = 1. The control parameters Ap and np, with p = k — m + 1, k — 
m + 2,..., k, are dehned by the asymptotic boundary condition 

- Gm{x) {x^oo) . (100) 

Consequently, the self-similar root approximants can always be constructed, even when 
the number of terms in the left and right asymptotic expansions are not equal to each 
other. 
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6 Self-Similar Exponential Approximants 


A different strategy is to be pursued when only a single-side asymptotic expansion, say 
at a; —0, is available. Then there are no boundary conditions determining control 
parameters. The latter are to be specihed in a different way appealing to convergence 
properties. A particular choice of control parameters results in a nice structure of nested 
exponentials [70,72,75]. Here we present a more rehned derivation of the exponential 
approximants and suggest some novel ways of constructing the cost functionals dehning 
the effective approximation time, that is, the control time. 

Looking at the fractal transform (71), it is easy to notice that the convergence prop¬ 
erties of the sequence {ids} improve if |a;|® —> 0. The latter can be realized when a; —0, 
i.e. in the same situation as for the asymptotic expansion (66). But there are in addition 
two other possibilities for |a;|® to tend to zero, when 

-|-cx) , |a;| < 1 

—cx) , |a;| > 1 

By taking these limits, we may extend the region of applicability of the function presented 
by the asymptotic expansion (66), valid only at a; —0, to the regions [0,1) and (1, cx)). 

To derive a self-similar approximant for the case (101), by employing the method of 
fractal transforms, we need to obtain a kind of an iterative equation, similar to Eq. (83). 
For this purpose, let us introduce a set of functions 



( 101 ) 


iJk{zn) = l + Zn {n = l,2,...,k) , 

Zn Zn{x^ , (^fc-|-l) = 1 , 

being iteratively connected with each other by means of the relation 

Z^{x) = ^ 

1 

Then the series (66) can be identically presented as 

gk{x) = i)k{zi) . 


( 102 ) 

(103) 


(104) 


The self-similar renormalization of ipkizn), accomplished by means of the method of fractal 
transforms, is 

'lp*k{Zn,s) = (^1- ^ TnZn^ ■ (105) 

Realizing the /c-step renormalization for the iterative relations (102) and (103), we hnd 

^fc(a;,s) = ^7fc«,s) , (106) 

where z* = z^{x, s), with n = 1, 2,..., /c, and 

Zn{x, s) = ^fc(4+i, s) . (107) 

^n—1 

According to the last identity in Eqs. (102), 


( 4 + 1 , 5 ) 


1 . 
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Irrespectively to what limit, either s +cx) or s —> —oo, is taken in the form (105), 
one gets the same result 

lim (l — - = expfr;^) . 

s^itoo \ s J 

Therefore, in what follows, we may write |s| — oo, keeping in mind any of the limits 
s —> ±cx). Let us dehne 

= lim . (108) 

|sj^oo 

And introduce the notation 

Cfi = Tfi , iy>fi = (Xfi (Xfi—i . (109) 

^n—l 

Taking the limit (108) in the iterative relations (106) and (107), we come to the self-similar 
exponential approximant 

gl{x) = exp {cix'^^ exp {c 2 x''^... exp {ckx''’")).., (HO) 

for short, called superexponential. 

Expression (110) contains the coefficients c„ that, as is seen from the notation (109), 
are proportional to the control time r^, which is not yet dehned. The simplest way would 
be to set Tn = 1/n, as in Eq. (53). It could also be possible to find from a fixed- 
point condition. However, the most general and rehned way is to determine the control 
time by minimizing a cost functional [75]. In optimal control theory, one constructs 
cost functionals by formulating the desired properties of the system. For our case, the 
procedure can be as follows. If it is recalled that the control time describes the minimal 
time necessary for reaching a fixed point at the n-th step of the calculational procedure, 
then uTn approximately corresponds to the total time r* required for reaching the hxed 
point. When ~ r*, this implies that Tn ~ r*/n. The time of reaching a hxed point 
depends on how far this point is. The shorter is the distance from the point, the faster 
is the way to it. The distance passed at the n-th step can be evaluated as vfTn, with vf 
being a characteristic velocity at this step. In this manner, we need to hnd a minimal 
time Tn that is close to t* /n and which corresponds to the fastest passage of the distance 
v^Tn. These requirements suggest to construct the fastest passage cost functional 

= (h- —) +C{Krnf , (111) 

in which the parameter C ^ 0 is included for generality. The value of ( can be chosen 
if some additional information on the system is available. In the absence of such an 
additional information, we set C = 1- 

Dehning the characteristic velocity v*, it is natural to associate it with a cascade 
velocity Vn(x) taken at the most dangerous value of x, where convergence is the worst 
and, respectively, the deviation v*Tn should be the largest. Thinking back to the fractal 
transform (71), we know that the sequence {F^} converges, under condition (101), if either 
|a;| < 1 or |a;| > 1. This means that the dangerous point is |a;| = 1. Therefore, we dehne 

v* = Vnix) (|a;| = 1) . (112) 
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With the cascade velocity 


Vn{x) = gn{x) - gn-l{x) = OnX"" , 


we have 


(vl)^ = ai 


So that the fastest passage cost functional (111) becomes 


^r = -y\ 
2 ^ 


r 


* \ 2 


r, 


n — ) + Calrl 
n 


(113) 

(114) 

(115) 


The control time is given by the minimization of the cost functional (115), i.e. from 
the conditions 

= 0i ^>0- (116) 


SXn 

The extremum condition leads to 




r 


r„, = 


n(l + Ca^; 

The found extremum is a minimum, since 

52 


(117) 




— 1 + C®n > 0 • 


What has been yet left undehned is the effective total time r*, which can be derived 
from the following reasoning. If the sought hxed point is reached in one step, this implies 
that Ti = 1. Applying this condition to formula (117) yields 

T* = 1 + (af (n = 1) . (118) 


Then the control time is 


r„= (119) 

n{l + Cal) 

The parameters c^, dehned in Eq. (109), are proportional to the control time r^, because 
of which they can be called the control parameters or, simply, controllers. With the 
control time (119), the controllers are 


a„(l + Cai) 

na„_i(l + C^n) 


( 120 ) 


In this way, the superexponential (110) is completely dehned. 

If the function g{x) was introduced as a scale-invariant form of the sought function 
/(x), then the self-similar exponential approximant for the latter is fk{.x) = ipQ{x)gl{x). 
Recall that the function g{x) has been assnmed to be hnite on the manifold X. In the 
case of a fnnction g{x) divergent at some point Xq G X, one should consider its inverse 
g~^{x), provided this is a hnite everywhere on X. In the example of deriving the control 
time (119) from the cost functional (115), it is snpposed that the function f\x) is sign 
dehnite so that the function g{x) is nonnegative. When it is known that f{x) changes 
its sign, this information has to be encompassed in the procedure. This can be done, 
for example, by factoring f{x) = (p{x)g{x), with g{x) being positive. Another possibility 
could be to incorporate information on the points of the sign change into the constructed 
cost functional. For describing oscillating fnnctions, it could be conceivable to deal with 
complex control times. 


23 



7 Examples 

Self-similar approximation theory has been applied to various physics problems. These 
applications can be found in the cited references [35-57,68-75]. Among recent works, 
we may mention the usage of this approach to barrier crossing processes [78,79], critical 
phenomena [80], and to the rupture of mechanical systems [81]. In the present section, 
we give several examples which illustrate some new possibilities of the approach. 

A. Amplitude of Elastic Scattering 

This example is interesting by demonstrating the use of superroots when the number 
of terms in the left asymptotic expansion is much larger than that in the right expansion. 

Consider the scattering of two particles of masses mi and m2, with momenta pi and 
P2 before collision and p'^ and P2 after it. The four-momenta are normalized on the mass 
shell so that pj = m]. The scattering amplitudes are usually presented as functions of 
the Mandelstam variables [82] which are 

s = {pi+ P 2 f = {Pi + P 2 f , t = {Pi- P'lf = {P 2 - P 2 f , 

u = {pi- p'^f = {p2 - p'l^ , s + t + u = 2 (mf + ml) . 

The amplitude of elastic scattering can be expressed, be means of perturbation theory, as 
an asymptotic expansion in powers of the coupling parameter g, 

00 

Tig,s,t) g + J2Tn{s,t)g^ (g ^ 0) , (121) 

n =2 

where Tn{s,t) ~ exp(ns) as s —cx). It is known that for any g there exists the Froissart 
upper bound given by the inequality \T{g,s,t)\ < |A(s,f)|, where a particular form of 
A{s,t) depends on whether the considered theory is local [82] or nonlocal [83]. Since the 
Froissart upper bound is valid for any g, including cx), let us assume that 

T{g,s,t) A{s,t) (^^cx). (122) 

Our aim is to construct a crossover formula between the left and right expansions (121) 
and (122), respectively. 

Following the general scheme, we, first, have to introduce the scale-invariant function 

f{g) = -T{g,s,t) , (123) 

9 

in which, for brevity, we do not write explicitly other variables, except g. Denoting for 
the function (123) the weak-coupling. 


fig) - Pkig) 

ig 

^0), 

(124) 

and the strong-coupling, 

fig) ^ ^i(^) 

ig - 

00 ), 

(125) 

asymptotic expansions, from Eqs. (121) and (122), 

we have 


V^fc(fi') — 1 + o-i9 + + ... 

+ dkg^ 

' , ^i(^) = bg~^ , 

(126) 
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where = Tn+i{s,t) and b = A{s,t). 

The crossover formula for the scattering amplitude 

Tkig) = g^hig) 

is obtained by constructing the self-similar root approximant (99) for 

vlM = + As'')”* , 


( 127 ) 


(128) 


where the control parameters and A^. are defined by the asymptotic boundary condition 
(100), which gives 

nfc = - ^ , Ak = ^ . (129) 

Thus, a self-similar root approximant for the scattering amplitude is 


r;(9) = 


!()'■ +9>Li(9)9‘I‘''‘ ' 


(130) 


Note that the expansion (p^-iig) here can also be converted to given by the 

superexponential (110). 

B. Summation of Numerical Series 

Poorly convergent or divergent numerical series can be summed by means of the su¬ 
perexponentials in the following way. Let us consider a series Soo = Y^^=o(^n, whose 
particular sums are 


A ^ XI 9 n ■ 


(131) 


Introduce the function 




(132) 


for which S'fc(l) = Sk- Construct the self-similar exponential approximant Sl{x), accord¬ 
ing to formula (110). Setting a; = 1 in S'^(a;), we get S]: = which is the sought 

self-similar approximant 


Sj = ao exp (ci exp (C2 ■ ■ ■ exp(ci))...) , 


(133) 


with the controllers c„ given by Eq. (120), where we set C = 1- If fh® exact value of S^o 
is known, one can compare the accuracy of the particular sums (131), characterized by 
the percentage error 

= (134) 


( 135 ) 


with the accuracy of the self-similar approximants (133), described by the error 


:l= -A- - 1 ■ 100 % . 


As an illustration, let us consider the sum (131), with the coefficients 


2n -I- 1 
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The sequence {Sk} converges to Soo = 7r/4. This convergence is rather slow, for example, 
the percentage errors (134) for the hrst hve terms are 

-15%, 10%, -7.8%, 6.3%, -5.3%, 

while the superexponential (133) gives the errors (135) for the hrst hve approximants as 

-8.8% , -3.6% , -1.9% , -1.6% , -1.6% , 


demonstrating a much faster convergence. 

C. Multiloop Feynman Integrals 

Employing the Feynman diagram techniques in quantum held theory or quantum 
statistical mechanics, one confronts with the so-called multiloop integrals. These can be 
calculated by means of perturbation theory [84,85] resulting in asymptotic series. The 
latter can be summed with the help of the superexponentials. 

Let us start the illustration with a simple one-loop integral 


where a is a positive parameter, D is space dimensionality. The exact value of the 
(136) is known to be 


I{a,D) 


r{a-D/2) 
(47r)'^/^r(a) ’ 


integral 

( 137 ) 


where r(-) is a gamma-function. A perturbative procedure for Eq. (136) can be dehned 
[84,85] by introducing 


/(a, D, e) 


1 

(27r)^ 


d^p 

(1 -F £p2)“A 


(138) 


and expanding the integrand in powers of e, which, after the integration term by term, 
results in a series 


I{a,D,e) = . 

n 


( 139 ) 


Since, as follows from Eqs. (136) and (138), I{a,D, 1) = I{a,D), the answer is obtained 
by setting e = 1. 

Accomplishing a partial self-similar exponentization of the series (139), one gets 


with 

g{e) = Jeans'" , 

n 

the coefficients a„ being 

_D + l _D{D + 2) 

^ “ “ 6 ^ ’ “ 240 ^ ’ 

_{D + 1)(3T)2 + 6T) - 4) _ D{D + 2){D‘^ + 2D-2) 

“ 240a3 ’ “ 240a4 


(140) 


( 141 ) 
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The partial sums of Eq. (141), after setting e = 1, become 


9k 


( 142 ) 


The corresponding superexponentials are 

gl = exp (ci exp (ca ... exp(cfc))) . 

Finally, for the integral (136), we hnd the self-similar approximants 


4 *(a,^) = Ta — 

(47ra)^/^ 


D{D + 2) 


9k-i r ’ 


(143) 


(144) 


with Qq = 1. 

Consider the case of a = 1 and D = 1, when /(1,1) = 1/2. The perturbation series 
(139) take the form 


1 1 1 / 3 25 , 


105 _3 1659 ^4 

1024 ^ 32768 


(145) 


The coefficients a„ in the sum (142) are 


_ 

ao = 1 , ai = - , 02 = - , ^ 3Q 


For the self-similar approximant (144), we have 


44/F (s 


(146) 


The errors of the perturbative expression (145) at e = 1 are 

-44% , -22% , -11% , -5.6% , -2.8% , 

which is to be compared with the errors of the self-similar approximants (146), 


-7% , -4. 


-0.77% , -0.14% , -0.08% , 


which are an order smaller. 

For the case a = 2, D = 3, one has 1(2, 3) = l/87r. The perturbative expression (139) 
reads 


1(2,3,e) ~ 


15 385 2 4725 o 228459 

1 + 777^ + TTX ^ ^ + TTTTXXX 


' ’ ’ ' (87r)3/2 V 16 512 8192 524288 

The coefficients a„ from Eq. (142) are 

_ 1 _ 5 _ 41 _ 13 

00 - 1 . o.--, a,-- a,- — , a,- — . 

The self-similar approximant (144) becomes 


(147) 


5(2,3) = 


(8ir)3/2 


exp 77 9t-i 


(148) 
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The direct expansion (147) yields the errors 


-80% , -61% , -46% , -35% , -26% , 


while those of the self-similar approximants (148) are lower: 


-30% , -26% , -16% , -13% , -12% . 


Now, let us turn to a H-dimensional three-loop Feynman integral 


J{D)^ 


d^Pi d^p2 d^ps 


jD 


D, 


{2nr^J (l + pf)(l + p2)(l+p|)[l + (p, + P2 + P3)2] • 


(149) 


Following the same procedure as in the calculation of the previous Feynman integrals, 
one dehnes J{D,e) and then set e = 1. For concreteness, let us take D = 2. Then the 
self-similar approximants for integral (149) are dehned as 


JU2) = 


2567r=> 


exp 


9 


9k-i 


(150) 


with gl having the form (143), which is obtained from of Eq. (142), where the coefficient 

Ojyi cir6 


Oo — 1 , 


05 


7 

13 

59 



do —- do — 

144 

768 ^ 

1 ^4 — 

2324 

15243 

07 = 

150379 

18432 ’ 

“ 86016 

393216 


373 
3^ ’ 


The accuracy of the approximants (150) again is much better than that of simple pertur¬ 
bative expressions. The first seven approximants demonstrate a fast increase of accuracy. 
The related errors, calculated by comparing the self-similar form (150) with the numerical 
value J(2) = 0.00424027, are 


-46% , -40% , -29% , -24% , -21% , -20% , -19% . 

This demonstrates a monotonic convergence, while the standard perturbation theory in 
powers of e wildly diverges. 

Thus, the self-similar approximants provide rather good approximations even for very 
bad, fastly divergent series derived by means of standard perturbation theory. 

In conclusion, we may mention that the self-similar approximation theory has been 
successfully applied not only to a number of physical problems [35-57,68-75,78-81] but 
also to other complex systems, such as hnancial markets [86-90]. Time series, related to 
hnancial, economic, biological, and social systems, are known to possess special fractal 
properties [91-95]. This is why these series can be naturally described by self-similar 
approximants, especially by those that explicitly display their self-similar structure as in 
self-similar roots and self-similar exponentials. 
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